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ABSTRACT 

About 54% of radio galaxies at z ^3 and 8% of radio galaxies at 2< z <3 show 
unusually strong Lya emission, compared with the general population of high red- 
shift (z >2) radio galaxies. These Lya-excess objects (LAEs) show Lya/Hell values 
consistent with or above standard photoionization model predictions. 

We reject with confidence several scenarios to explain the unusual strength of 
Lya in these objects: shocks, low nebular metallicities, high gas densities and absorp- 
tion/scattering effects. We show that the most successful explanation is the presence 
of a young stellar population which provides the extra supply of ionizing photons 
required to explain the Lya excess in at least the most extreme LAEs (probably in 
all of them). This interpretation is strongly supported by the tentative trend found 
by other authors for z ^3 radio galaxies to show lower UV-rest frame polarization 
levels, or the dramatic increase on the detection rate at submm wavelengths of z >2.5 
radio galaxies. The enhanced star formation activity in LAEs could be a consequence 
of a recent merger which has triggered both the star formation and the AGN/radio 
activities. 

The measurement of unusually high Lya ratios in the extended gas of some high 
redshift radio galaxies suggests that star formation activity occurs in spatial scales of 
tens of kpc. 

We argue that, although the fraction of LAEs may be incompletely determined, 
both at 2< z < 3 and at z ^ 3, the much larger fraction of LAEs found at z ^3 is 
a genuine redshift evolution and not due to selection effects. Therefore, our results 
suggest that the radio galaxy phenomenon is more often associated with a massive 
starburst at z >3 than at z <3. 

Key words: cosmology: observations, early Universe; galaxies: active; galaxies: evo- 
lution 



1 INTRODUCTION 

High redshift powerful radio galaxies (HzRG, z >2) are char- 
acterized by an emission line spectrum which is rich in emis- 
sion lines of different ionization species. In the optical (UV 
rest frame), Lya is usually the strongest emission line, fol- 
lowed by CIVA1550, HeIIA1640, CIII]A1909 and, less fre- 
quently, NVA1240 (CIV, Hell, CIII] and NV hereafter). 



* Unpublished results are presented in this paper for the radio 
galaxy 1338-1941, based on observations carried out at the Euro- 
pean Southern Observatory, Paranal (Chile) for the ESO project 
69.B-0078(B). 



De Breuck et al. (2000a, DBOOa hereafter) observed that 
Lya is unusually strong in some HzRG, most of which are 
at z >3. They also found that the Lya equivalent width in- 
creases at z > 3 and the ratios of Lya to CIV, CIII] and Hell 
are roughly twice the value found for lower redshift objects 
(2< z < 3). The authors interpret these results as spectro- 
scopic evidence for the youth of radio galaxies at z >3, which 
are still surrounded by large halos of primordial hydrogen 
from which the galaxy is formed. 

In this paper we present a more thorough study of 
DBOOa] results. We start by quantifying (§3) the variation 
of the Lya/Hell and Lya/CIV ratios and Lya luminosity 
with redshift for radio galaxies at 1.8^ z ^4.4 and inter- 
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pret the results in the light of appropriate photoionization 
models (described in §3) and/or scenarios (§5). The results 
and the cosmological consequences will be discussed in §6. 
Summary and conclusions will be presented in §7. 

A Aa = 0.73, D. m = 0.27 and H = 71 km s" 1 Mpc" 1 
cosmology is adopted in this paper. 



2 OBJECT SAMPLE. 

The object sample has been extracted from the data set of 
DBOOa and De Breuck et al. (2001, DB01 hereafter), which 
includes the largest compilation of line fluxes for z >2 radio 
galaxies available. We will not consider lower z objects, for 
which Lya is not observable. Our final sample consists of 13 
radio galaxies at z ^3 and 48 radio galaxies at 1.8^ z <3. 

Table 1 shows the sample with some relevant informa- 
tion. 

The original data set of De Breuck et al. includes radio 
galaxies from nine samples designed to find HzRG. As the 
authors explain, the surveys can be divided into two types 
(see DBOOa for all appropriate references of each survey): (i) 
several large flux density-limited surveys such as the 3CR, 
MRC, BRL sample, and 6C, 7C and 8C surveys; (ii) several 
"filtered" surveys, which were designed to find the highest 
redshift objects. In the latter case, the radio spectral in- 
dex flux is most often used (ultra steep spectrum sources), 
sometimes in combination with an angular size upper limit. 
Alternatively, the filter consists of a flux density interval 
centered around the peak in the source counts around 1 Jy. 

The total sample contains 176 radio galaxies with z in 
the range 0.05 to 5.2. It contains a similar number of sources 
at z < 2 and at z >2. For the present study, we have elimi- 
nated ~110 objects for the following reasons: 

• 4 broad emission line objects (FWHM>4000 km 
s _1 and/or additio nal qua sar-like properties). These are 
2313+ 4053 (2 =2.99. lDB0ih , 0744+464 (z =2. 93, M cCarthy 
Il99ll ) 2025 -218 (z =2.63, Villar-Martm et al. Il999l . Larkin 
et al. |2000| ) and 0349-211 (2 =2.33, McCarthy et al. ll99ll ). 
For these objects, the UV lines have a contribution from 
the broad line region, where the physics is complex and not 
interesting for the current study. 

• 2 objects with broad absorption l ine sys tems: 1908+722 
(2 =3.53) and 1115+5016 (2 =2.54) (|PB0ll) . 

• 91 objects with no Lya measurement. Except for 3 ob- 
jects for which the reason is not known, in all other cases 
Lya is outside the observed spectral range. 

• 23 objects with no measurement of CIV and Hell. The 
reasons are varied: for 7 objects, the CIV and Hell lines 
were outside the useful spectral range, or only Lya narrow 
band imaging fluxes are available. For 6 objects, the lines are 
reported to be detected, but no fluxes are provided, except 
for Lya. For 7 objects, the lines where within the observed 
spectral range, but not detected. In 2 cases Hell was outside 
the observed spectral range and CIV was inside, but not 
detected. For 1 object, the reason for the non-availability of 
the CIV and H ell fluxes is unknown. 

Using iDBOll original spectra, we have been able to mea- 
sure useful upper limits (in some cases the spectra are too 
noisy or severely affected by sky residuals) for the CIV 
and/or Hell lines of several of these objects. The lower limits 
for the Lya ratios are shown in Table 1. 



There are 10 objects with no measurement of CIV (2) 
or Hell (8). In all cases, except one, both lines where within 
the observed spectral range, but not detected. Upper limits 
could be measured for two objects using the original DB01 
spectra. The lower limits for the Lya ratios are given in 
Table 1. 

The following objects were part of De Breuck et al. sam- 
ple, but the line fluxes of one or more of the lines of interest 
to us were not provided. For this reason, we present measure- 
ments based on our own long slit Keck LRISp spectroscopic 
data or results published by other authors. 

• 4C40.36 at 2 =2.27 (nuclear aperture, Humph rey |2005l . 
Humphrey et al. 2006, in prep; Vernet et al. l200ll ) ) 

• 4C48.48 at 2 =2.34 (nuclear aperture, Humph rey 120051 . 
Humphrey et al. 2006, in prep; Vernet et al. l200ll ) 

• 731+438 at 2 =2.43 (nuclear aperture, H umph rey 
120051 . Humphrey et al. 2006, in prep; Vernet et al. 120011 ) 

• 0902+34 at 2 =3.40 (Martm-Mirones et al. 1 19951 ) 

• 1243+036 at 2 =3.56 (nuclear aperture, H umph rey 
l2005l . Humphrey et al. 2006, in pre p; Ver net et al. l200ll ) 

• 4C41.17 at 2 =3.79 (Dey et al. Il997l ) 

• 1338-1941 at 2 =4.11 (De Breuck et al. in prep. ESO 
project 69.B-0078(B).) 

More accurate line measurements than De Breuck et al. 
original data were used for: 

• 2334+1545 at 2 =2.48. Hell is not detected in this ob- 
ject. Because of t he unu sually narrow CIV line (<500 km 
s~ ) measured bv lDBOll compared with Lya (1900 km s -1 ) 
and CIII] (1600 km s" 1 ) lines, we checked this spectrum 
again. We cannot confirm that the line is detected. We con- 
clude that Lya/CIV>5 and Lya/HeII>5. 

• 2104-242 at 2 =2.49 ( based on VLT FORS long slit 
spectroscopy, Overzier et al. l200ll) 

• 0140+3253 at 2=4.41. According to IDBOll. CIV is de- 
tected and Lya/CIV=43±10. Because of the strong sky 
residuals in the CIV region, we have checked the spectrum 
again. We cannot confirm the detection of the line. We have 
measured upper limits for both CIV and Hell. We obtain 
Lya/CIV>30 and Lya/HeII>10. 

• Upper limits for CIV and/or Hell were measured 
for 0231+3600 (2=3.08), 0040+3857 (2=2.61), 2355-002 
(2 =2.59) and 2254+185 (2=2.15) using the original spectra 
of lDBOll . 



3 MODELS 

Photoionization models were computed using the multipur- 
pose code Mapp ings Ic (Binette, Dopita & Tuohy [l985l . Fer- 
ruit et al. Il997l ). The following assumptions will be applied 
unless stated otherwise: the ionizing continuum is a power 
law (PL) of index a=-1.5 (F v ~ v a ) with a cut off en- 
ergy of 5xl0 4 eV. The gas is characterized by a density 
n =100 cm -3 and solar abundances. The clouds are radia- 
tion bounded and the geometry is plane parallel. The den- 
sity behaviour within the clouds is isobaric. Dust reddening 
effects have been ignored in the modeling. These are the 
standard models which, in general, reproduce most of the 
UV and optical line r atios of HzRG rather succe ssfully (e.g. 
Robi nson et al. Il987l . Villar-Martm et al. Il997l . Humphrey 
120051) . 
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Object 


z 


Lya/HcII 


Lya/CIV 


P(%) 


S850/xm (mjy) 


0140+3253 


4.41 


>10 


>30 


— 


3.3±1.5 


1338-1941 


4.11 


18+4 


11+2 


5.0 


7±1 


4C60.07 


3.79 




3.7 




17±1 


4C41.17 


3.79 


27±2 


11.1+0.7 


<2.4 


12.1±0.9 


1911+6342 


3.59 


18+2 


18+2 


- 


<11.9 


2141+192 


3.59 


17.8 


10.6 


- 


<5 


1243+036 


3.57 


34±3 


14+1 


11±4 


<6.5 


0121+1320 


3.52 


7.1+2 


8+2 


7±2 


7.5±1.0 


0205+2242 


3.51 


26+4 


18+2 


- 


<5.2 


0902+34 


3.40 


9.1 


6.3 


1.5 


- 


1123+314 


3.22 


2.6+0.4 


1.8+0.3 


- 


5±1 


1112-2948 


3.09 


2.4+0.3 


4.1+0.7 


- 


5.9±1.6 


0231+3600 


3.08 


>5 


>5 


- 


6±2 


0747+365 


2.99 


3.5+0.7 


4.6+0.9 




5±1 


0943-242 


2.92 


7.4 


5.2 


6.6±0.9 


— 


4C24.28 


2.88 




4.3 


— 


2.4±1.2 


0857+036 


2.81 


3.7 


2.6 


— 


— 


0417-181 


2.77 


6.0 


2.5 


— 


— 


1019+0534 


2.77 


1.0 


0.8 


— 


<3.8 


1338+3532 


2.77 


5.9+0.8 


14+3 


— 


— 


0920-071 


2.76 


15+2 


12+2 


— 


— 


1545-234 


2.76 


3.7 


6.1 






2202+128 


2.71 


3.2 


4.1 


— 


— 


1357+007 


2.67 


— 


5.4 


— 


— 


0040+3857 


2.61 


>5 


4+2 


— 


— 


2355-002 


2.59 


>3 


4.0+1 


— 


— 


0529-549 


2.58 


12.3 


18.5 


— 


— 


0828+193 


2.57 


7.0 


7.0 


10±1 




1755-6916 


2.55 


23+4 


15+2 


— 


— 


1436+157 


2.54 


7.0 


2.5 


— 


— 


1558-003 


2.53 


8.8 


5.5 


— 


— 


0303+3733 


2.51 


21+4 


7+1 


— 


— 


1650+0955 


2.51 


8+1 


7+1 


— 


— 


2104-242 


2.49 


10+1 


13+2 


— 


— 


2334+1545 


2.48 


>5 


>5 


— 


— 


4C23.56 


2.48 


5.3 


3.8 


15±2 


<5.5 


2308+0336 


2.46 


7.5 


4.7 


— 


<3.2 


0731+438 


2.43 


14+1 


8.6+0.7 


<2.4 


— 


1033-1339 


2.43 


12+2 


4.3+0.7 


— 


— 


0748+134 


2.42 


4.2 


3.5 


— 


— 


1410-001 


2.36 


11.1 


8.3 


12±3 


— 


1707+105 


2.35 


4.8 




— 


— 


0128-264 


2.35 




7.2 


— 


— 


1428-1502 


2.35 




9+2 






4C48.48 


2.34 


■ 


8.2+0.8 


8.4±1.5 


5±1 


0211-122 


2.34 


1.0 


1.8 


19dbl 




1251+1104 


2.32 


7.7 


7.7 


— 


— 


4C40.36 


2.27 


11.6+0.8 


7.3+0.7 


7±1 


<4.1 


1113-178 


2.24 


9.1 


3.8 


— 


— 


0200+015 


2.23 


5.4 


4.1 


— 


— 


1138-262 


2.16 


10.7 


17.4 




13+3 


0355-037 


2.15 


3.0 


4.1 






2254+185 


2.15 


10+5 


>4.0 


<2.6 




2036+0256 


2.13 


10+3 


11+5 






1102-1651 


2.11 


2.1+0.5 


2.7+0.7 






1802+645 


2.11 




8.9 






1740+664 


2.10 




13.7 






0448+091 


2.04 


8.7 


10.1 






1805+633 


1.84 




6.8 






3C256 


1.82 


9.9 


10.4 


12±2 




3C294 


1.79 


10 


10 




<2.3 



Table 1. Sample of radio galaxies used in this study in order of decreasing redshift and with some relevant information. LAEs are 
highlighted in bold, z ^3 a nd z < 3 radio galaxies are separated by a horizontal line. See §2 f or ref erences of the optical data. P(%) 
values are from Vernet et al. l200ll. D ev et al. Il997l . De Breuck et al. (in prep.) and Jannuzi et al. 1 19951 . Submm data are from Archibald 
et al. l200ll and Reuland et al. Tioo4 No available measurements are indicated with Errors for P(%) and Ssso are provided when 
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All the stellar SEDs (Spectral Energy Distributions) 
us ed in t his paper were built using Starburst99 (Leitherer et 
al. ll999h . a web based software and data package designed 
to model spectrophotometric and related properties of star- 
forming galaxies. We have assumed an initial mass function 
(IMF) with spectral indexes 1.3 and 2.3 for the mass inter- 
vals 0.1 to 0.5 M Q and 0.5 to 120 M Q respectively. Stellar 
abundances are always considered to be solar. Both contin- 
uous and instantaneous star forming histories were consid- 
ered. 



4 VARIATION OF Lya/Hell, Lya/CIV AND 
Lya LUMINOSITY WITH REDSHIFT 

In this section we investigate how the Lya/Hell and 
Lya/CIV ratios and the Lya luminosity change with red- 
shift. Both ratios have been measured in many radio galaxies 
at z >2, since the lines are observable in the optical window. 

Three diagnostic diagrams involving the Lya, Hell and 
CIV lines are shown in Fig.l. The solid line represents stan- 
dard photoionization models (see §3). The ionization param- 
eter varies within the range [0.0046,1]. 

The Lya/Hell ratio predicted by the models is in the 
range ~15-20. Due to its strong sensitivity to U, the pre- 
dicted Lya/CIV values cover a much larger range. The min- 
imum predicted value is ~9. Because Lya/Hell shows little 
dependence on U compared with Lya/CIV, we will focus 
most of our discussion on the first ratio, with reference to 
other ratios when useful. 

All objects in Table 1 with measurements available for 
the three line ratios (including lower and upper limits) are 
plotted in Fig. 1. We have marked with coloured symbols 
those radio galaxies with Lya/Hell values consistent with 
or above the standard case B value (0140+3253 at z =4.41 
is the only exception, see below). Given that errors are not 
available for many radio galaxies in the sample, we con- 
sider an object to be consistent with the model predictions 
when Lya/Hell is at least 90% of the minimum values pre- 
dicted by the models, i.e. Lya/HeII^14. Because of the very 
high values compared with typical values of most HzRG, we 
will call these radio galaxies Lya excess objects (LAEs here- 
after). Black symbols correspond to objects with Lya/Hell 
ratios below the case B recombination value. 

The diagrams show a very clear segregation of the most 
distant radio galaxies in the sense that most z ^3 objects 
have the highest Lya/Hell values. The Lya/CIV ratios of 
these objects are also among the largest values. 

Of the 10 sources at Js3 for which Lya/Hell has been 
measured, 6 have Lya/Helljsl4 and are therefore LAEs. 4 
are not LAEs. It is not possible to classify the remaining 
three z ^3 radio galaxies due to the unavailability of the 
Lya/Hell value. 

One more z ^3 radio galaxy can be added to the group 
of LAEs: 0140+3253 (z =4.41), which shows Lya/HeII>10 
and Lya/CIV>30. The large Lya/CIV ratio does not nec- 
essarily imply that the object is a LAE, since this ratio can 



1 U = 47r ^ ne , where Q is the photon ionizing luminosity, r is 
the distance between the cloud and the ionizing source, n is the 
particle gas density and c is the speed of light. 



be reproduced by very low or high U values (see Fig. 1). 
The real value of Lya/Hell would be necessary. However, 
the fact that this object presents the largest Lya/CIV value 
of the whole sample suggests that this is also a LAE in the 
sense that the Lya emission is unusually enhanced. This is 
also suggested by its high Lya luminosity (1.35xl0 44 erg 
s _1 ), which is more characteristic of LAEs than non LAEs 
(see below). Therefore, we classify 0140+3253 as a LAE. 

We conclude that at least 7 out of 13 (54%) z ^3 radio 
galaxies are LAEs. On the other hand, only 4 out of 48 (8%) 
2< z <3 radio galaxies can be classified as such (see Table 
10. 

We show in Table 2 the median and average (1st num- 
ber in brackets) values of the Lya ratios for the two redshift 
samples together with other interesting information (upper 
and lower limits have been excluded in the calculations). 
Lya/Hell and Lya/CIV are ~2.2 and ~1.5 times larger 
in z ^3 radio galaxies, which also tend to have noticeably 
larger Lya (~3 times) luminosities compared with the low z 
sample. The CIV/Hell ratio is similar in both groups. The 
bottom line shows the probabilities P that the two groups 
are drawn from different parent populations according to the 
Kolmogorov-Smirnov test. These values confirm that z 5=3 
radio galaxies show larger Lya ratios and higher Lya lumi- 
nosities. Interestingly, there is no distinction in radio size 
between the two samples. 

We have repeated the calculations from Table 2, this 
time separating the sample into LAE and non-LAE sources, 
rather than into z <3 and z 5s3 sources. The results are 
shown in Table 3. LAEs tend to be at higher redshift with 
a median z value z m ed~3.57, while z rne d=2.48 for objects 
with no Lya excess. They have Lya luminosities ~3.4 times 
higher and Lya/Hell and Lya/CIV ratios 2.7 and 2.1 times 
higher respectively compared with non-LAEs. Interestingly, 
also the CIV/Hell ratio is 1.6 times higher. There is also a 
clear difference in radio size. The bottom line shows the P 
values resulting from the Kolmogorov-Smirnov test. 

Interestingly, there are 4 LAEs (out of 11, i.e. 36%) 
which, even after taking into account the errors, show a 
Lya/Hell well above the maximum value (20) predicted 
by the models. These objects are: 1338-194 1 (Lyq /HeII~30 
once corrected for absorption, Wilman et al. l2004T ). 4C41.17, 
1243+036, 0205+2242, all at z >3. No or very weak absorp- 
tion has been found for the LAEs 4C41.17, 1243+036 (van 
Ojik et al. 1997) and 1755-6916 (Wilman et al. 2004). This 
information is not available for the rest of the LAEs. 



4.1 Measurements of spatially extended Lya 
ratios 

The results above are based on spatially integrated spectra, 
centered on the high surface brightness regions (although 
not necessarily covering the whole extension of the emission 
line regions). There are two z ^3 objects in the literature 
for which the spatial variation of the Lya ratios has been 

2 Considering a power law of index a=-1.0 (Villar-Martin, 
Tadhunter & Clark, 1997) would not change our conclusions. 
Lya/Hell would be in the range ~10-14, implying that still the 
vast majority of 2< z <3 objects lie below the model predictions, 
while most z ^3 radio galaxies are LAEs. 
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z range 


Lya 
Hell 


Lya 
CIV 


CIV 
Hell 


L(Lya) 
10 44 erg s — 1 


LAS 


2< z < 3 


8.00 (8.32,[4.75]) 


6.90 (7.24,[4.21]) 


1.20 (1.29,[0.65]) 


0.44 (0.74,[0.79]) 


52 (66.2, [60.6]) 


z >3 


17.50 (16.01, [10.90]) 


10.60 (10.07, [5.41]) 


1.44 (1.51, [0.59]) 


1.35 (1.34,(1.03]) 


43.0 (78.0,(74.7]) 


P 


99% 


96.5% 


71% 


93% 


1% 



Table 2. The median, average (first number in brackets) and standard deviation (2nd number in brackets) of the main UV line ratios 
and the Lya luminosities in the 2< z <3 and z ^3 samples are shown. The Lya/Hell and Lya/CIV ratios are ~2.2 and 1.5 times higher 
in the z ^3 sample. On the other hand, Lya is ~3 times more luminous in this sample. No distinction is apparent on radio size between 
the two samples. The bottom line shows the probability that the two samples are drawn from different parent populations according to 
the Kolmogorov-Smirnov test. (Radio largest angular size values (LAS, in Tables 2 and 3) taken from DB01, DBOOa and De Breuck et 
al. 2000b) 
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Lya/Hell Lya/CIV 

Power law, index -1.5, n=100, solar abundances, U sequence 



Lya/CIV 



Figure 1. Diagnostic diagrams involving the Lya, Hell and CIV emission lines (see electronic manuscript for colour version of the 
figures). The solid black line represents the standard U sequence of photoionization models (see text), with the U value indicated for 
some points along the sequence. The dashed arrows show the sense of increasing U. The Lya/Hell values predicted by the models are in 
the range ~15-20. Objects with Lya/Hell ratios consistent with or above such values (named LAEs throughout this work) are indicated 
in colour and their names shown in the adjacent labels. Error bars for these objects are shown when available. Radio galaxies at z <3 are 
shown with dashed errorbars and/or small sized symbols and radio galaxies at z JS3 are shown with solid errorbars and/or large sized 
symbols. Solid arrows indicate lower and upper limits. The Lya ratios have not been corrected for absorption effects, except for 1338-192 
for which both the absorbed (red solid circle) and corrected (red hollow circle) values are shown. The z JS3 sample tend to have larger 
Lya/HcII and Lya/CIV ratios. Lya/Hell is above the model predictions for several objects. 



investigated in the direction of the radio structures (both of 
them are in the sample above): 1243+036 (green solid trian- 
gle in the figures) at z =3.36 (Humphrey I2005I . Humphrey 
et al. 2006, in prep.) and 2141+192 (blue open square in the 
figures) at z =3.59 (Maxfield et al. l2002h . For 1243+036 the 
Lya/Hell ratio has values in the range 26+5 up to >40 in 
the outer parts (>7" from the nuclear region) of the object. 
Lya/HeII=33±l in the nuclear region. On the other hand, 
Lya/CIV presents a minimum in the nuclear region (14+1) 
and a maximum of >40 (this is a lower limit) in the outer 
parts of the object (>7"). The spatially integrated values 
are ~34±1 and 14+1 respectively. 

The sec ond object is 2141+192 at z =3.59 (Max- 
field et al. 2002). The integrated spectrum gives ~18 
and 11 for Lya/H e ll an d Lya/CIV respectively. Fig. 7 
in iMaxfield etafl (|2002h shows a region at ~1-1.5" 
from the continuum centroid, where Lya/HeII>40 and 
Lya/CIV~15-20. At ~4-4.5", Lya/HeII>30. 



iHumphrevI (120051 ) (Humphrey et al. 2006, in prep.) 
studied the spatial variation of the Lya ratios for a sample of 
10 radio galaxies (all included in our study) in the redshift 
range 2< z <3. Five of these objects (1558-003, 4C40.36, 
B3 0731+438, 4C48.48 and 4C-00.54) show Lya/Hell in the 
range 14-30 at some spatial positions in the extended gas. 
These large Lya/Hell values are measured in all cases out- 
side the nuclear region. Among these five objects, only B3 
0731+438 is classified as a LAE from its spatially integrated 
spectrum (blue hollow small square in the figures). 



5 SCENARIOS 

In this section, we will investigate a) why Lya is unusu- 
ally strong in LAEs relative to other emission lines and the 
continuum, compared with the majority of HzRG; 6) why 
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Object class 


z 


Lya 
Hell 


Lya 
CIV 


CIV 

Hell 


L(Lya) 
10 44 erg s" 1 


LAS 
kpc 


non-LAEs 


2.48 (2.53,[0.42]) 


7.25 (6.94,[3.30]) 


5.45 (6.70, [4.05]) 


1.00 (1.20,[0.60]) 


0.43 (0.77, [0.84]) 


85.2 (53.5, [74.7]) 


LAEs 


3.57 (3.35,[0.68]) 


19.50 (21.38J6.20]) 


11.50 (12.52, [3.69]) 


1.63 (1.80,[0.62]) 


1.45 (1.55,[1.05]) 


18.5 (38.8, [42.6]) 


P 


98% 




99.9% 


99.5% 


94% 


90% 



Table 3. As Table 2, but the objects have been classified as LAEs (radio galaxies with Lya excess) and non-LAEs. The LAEs tend 
to be at higher redshift. Lya is ~3.4 times more luminous in these objects and the Lya/Hell and Lya/CIV ratios are 2.7 and 2.1 
times higher than in non-LAEs. The radio size is ~4.6 times smaller for LAEs. The bottom lines shows the P values resulting from the 
Kolmogorov-Smirnov test. 



LAEs show Lya/Hell values which are often well above the 
standard model predictions (see §4). 

We are looking for a mechanism that can enhance the 
Lya emission over other emission lines and the continuum. 
We have considered two possibilities: different properties of 
the ionized nebulae (heavy element abundances, gas den- 
sity, absorption and scattering properties) and alternative 
ionization mechanisms (shocks, stellar photoionization vs. 
AGN photoionization). 

The diagrams in Fig. 1 show a large scatter of the mea- 
sured Lya ratios, which cannot be explained by the standard 
models. Most radio galaxies show too faint Lya relative to 
CIV and Hell. The most natural explanation is that the Lya 
photons are absorbed by neutral gas and dust. Therefore, 
the most realistic way to model the line ratios of the whole 
sample should take absorption effects into account. This in- 
formation is available only for a few objects (Wilman et al. 
2004). Van Ojik et al. (1997) found clear signs of absorption 
in 11 out of 18 radio galaxies at z >2, but they do not pro- 
vide the information necessary to correct the Lya fluxes in 
the absorbed objects. 

Our expectation is t hat LAE s should show stronger ab- 
sorption than non LAEs. iDBOOal found that Lya absorption 
occurs most frequently in so urces with small radio sources 
(see also van Ojik et al. Il997l ) and at higher redshifts. Since 
most LAEs have higher z and smaller radio sizes than non 
LAEs (see Table 3), it is likely that they are also more ab- 
sorbed. In this case, they are likely to be more overluminous 
in Lya than assumed in this work. 

On the other hand, the goal of this study is to explain 
the Lya excess in the LAEs. For this and the above reason, 
absorption will be ignored in our models, but only discussed 
qualitatively when necessary. It is, therefore, not to be ex- 
pected that the models presented in this section should re- 
produce the line ratios of the whole sample, but only those 
of the LAEs. 



5.1 Lower metallicities. 

We investigate next whether lower nebular metallicities in 
LAEs can explain the observational results. 

A decrease of the heavy element abundances will en- 
hance the nebular electron temperature, which will favour 
the collisional excitation of Lya (and heavy element lines 
such as CIVA1550), while this process will not affect 
HeIIA1640 or the continuum significantly. 

We have computed a set of photoionization models to 
investigate the impact of the metallicity variation on the 
Lya ratios. The density has been fixed to n=100 cm -3 . We 



show in Fig. 2 the same diagnostic diagrams as before, with 
two model sequences for U =0.05 (dashed black line) and 
U =0.015 (red solid line) respectively. A range of U values 
is considered, since the ionization level of the gas is known 
to vary from object to o bject (i.e., Robin son et al. Il987l . 
Villar-Martm et al. Il997l . Humphrey 120051 ). The metallic- 
ity range covered by the two sequences is 0.05 to 2 Zq. U 
values outside the quoted range would produce a) similar 
results (U—0.1 would produce similar results to (7=0.05, al- 
though implying even lower metallicities for the objects) or 
b) discrepant line ratios. For example, very low ionization 
models (U <0.005) produce too faint CIV emission so that 
CIV/HeII<0.4. On the contrary, very high ionization mod- 
els (U >0.3) which can explain Lya/HeII>25 predict too 
strong CIV, with CIV/HeII>4. 

As the metallicity decreases towards well below solar 
values, the UV line ratios (in particular Lya/Hell) change 
dramatically. For Z >50% Zq this ratio is rather constant, 
since the nebula is too cold to produce any noticeable colli- 
sional excitation of Lya. As the metallicity decreases further, 
the high electron temperature makes this mechanism very 
efficient and the Lya/Hell ratio increases by a factor of 2 in 
the metallicity range 5-50 % Zq. The models also succeed 
to reproduce the high CIV/Hell ratios for some LAEs. 

A u-turn in the model sequences is observed for the 
ratios involving CIV due to the relatively less efficient exci- 
tation of the line (compared with Lya and Hell) as the ions 
become more and more scarce as the metallicity decreases. 
Interestingly, metallicities as low as 5% still produce strong 
metal (CIV) lines relative to Hell, with CIV/HeII>l if the 
gas is highly ionized (U ~0. 05-0.1). 

5 LAEs are consistent with solar abundances, choos- 
ing the appropriate U value and taking errors into account 
(lower metallicities are not discarded, though, due to the U- 
Z degeneracy in this area of the diagrams). These objects 
are: 1911+6342 , 2141+192 and 1338-1941 at z >3 and 0920- 
071, 0371+438 at z <3. Once corrected for Lya absorption, 
one of them (1338-1941, red circles) shifts towards the low 
metallicity models (10% Zq). I.e. if absorption was taken 
into account, it is possible that more objects would require 
well below solar abundances. 

Including 1338-1941, 6 LAEs have Z <25% Zq accord- 
ing to their UV line ratios. Interestingly, very low abun- 
dances Z <10% Zq are inferred for 4 of these objects: 1338- 
1941, 1243+036, 4C41.17, 1338-1941 (once corrected for ab- 
sorption) and 0205+2242, all at z >3. 

According to these results, LAEs are, on average, char- 
acterized by lower abundances than radio galaxies with no 
Lya excess. Abundances of less than 25% Zq are implied in 
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Figure 2. Metallicity effects. Same diagnostic diagrams as in Fig.l. The black and red lines represent sequences of photoionization 
models for two different U values. The metallicity varies in the range 0.05 to 2 Zq in both cases. The black arrows indicate the sense 
of increasing metallicity. Metallicity effects are very successful at reproducing the position of LAEs in these three diagnostic diagrams. 
Notice that at least 7 out of 13 objects require abundances <25% Zq. However, the models fail to reproduce the strength of the NV 
emission in at least the most extreme LAEs (see text). 



~55% of LAEs. Since Lya absorption has been taken into 
account for only one object, the metallicities could be even 
lower. 

The NVA1240 line provides a strong test for the validity 
of the low metallicity scenario. According to the models, this 
line should be very faint in the most metal poor objects. Its 
flux has been measured in three of them: 1243+036, 4C41.17 
and 1338-1992. We show in Table 4 the predictions of the 
NV/Hell and NV/CIV ratios for the models that best re- 
produce the Lya ratios of these three objects, together with 
the measured values. It is clear that in all three cases the 
NV emission is observed to be much (2 to 4 times) stronger 
than expected from the model predictions. Nitrogen over- 
abundance has been frequently found in different types of 
objects. This is usually attributed to secondary production 
of nitrogen. However, this mechanism cannot work at the 
very low metall icities implied by the models (Vila-Costas & 
Edmunds Il993h . Therefore, these LAEs are unlikely to have 
low metallicities. 

We therefore conclude that low metallicities do not pro- 
vide a satisfactory explanation for at least the most extreme 
LAEs. 



5.2 Density enhancement. 

An efficient way to enhance Lya over other emission lines 
and the continuum is increasing the gas density. At low den- 
sities the cooling of the nebula is due mostly to some forbid- 
den lines (e.g. [OIII]A4959,5007, [OIL1A3727, etc), via radia- 
tive de-excitation of their upper energy levels in the original 
ion. As the density increases, the critical density will be 
reached for some of these transitions, so that de-excitation 
becomes collisional, rather than radiative and radiative cool- 
ing becomes much less efficient. The electron temperature 



rises as a consequence and Lya will have an increasing con- 
tribution of collisional excitation in the partially ionized 
zone, which does not affect the HeIIA1640 recombination 
line. 

In addition, when densities approach values ~10 5 cm" 3 , 
the fraction of energy emitted as Lya photons versus two- 
photon continuum emission due to the decay of the 2 2 S 
level of H tends toward un ity, ra ther than 0.67 at the low 
density limit (Binette et al. [l993). 

The dramatic increase of Lya/Hell with density can 
be clearly seen in Fig. 3, where the red-solid and black- 
dashed lines are sequences of photoionization models for two 
different U values (solid red line: £7=0.05; black dashed line: 
[7=0.7), with n varying along each sequence in the range 
10 1 " 7 cm -3 and 10 1-6 cm -3 respectively. CIV collisional 
excitation is also more efficient so that CIV/Hell increases 
dramatically with density, while Lya/CIV decreases. 

Fig. 3 shows that density effects (coupled with a vari- 
ation in U) can explain the location of all LAEs in the di- 
agnostic diagrams (intermediate U values would lie between 
the two sequences). It is important to note that high densi- 
ties not only help to reproduce large Lya/Hell ratios, but 
also the trend for LAEs to show larger CIV/Hell, well above 
standard model predictions in some cases (see Fig.l). There 
are however arguments that make this explanation fail. 

The implied densities for some LAEs (e.g. 1243+036 
and 4C41.17) are ~10 5 cm -3 , i.e. very high compared with 
the densities measured or inferred for the extended emission 
line regions of radio galaxies at all redshifts (< several to 
few hundred cm -3 ). On the other hand, such high densities 
exist in the nucl ear narrow line region (e.g. De Robertis & 
Osterbrock ll984l ). 

Is it possible that the emission line spectrum of LAEs 
is dominated by the high density nuclear gas? This could be 
the case for one of them, 1338-1941 (z =4.11, Lya/HeII~30, 
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Object 


Lya 

Hell obs 


Lya 
CIV obs 


NV 

Hell obs 


NV 

CIV obs 


(U, Z) mo d 


Lya 

Hell mod 


Lya 

CIV mod 


NV 

Hell mod 


NV 

CIV mod 


1243+036° 


34+3 


14+1 


0.9 + 0.1 


0.37+0.06 


(0.05, 5%Z Q ) 


32.3 


14.8 


0.26 


0.12 


4C41.17 6 


27+2 


11.1+0.7 


0.71+0.08 


0.29+0.03 


(0.05,10%Z o ) 


27.5 


10.0 


0.37 


0.13 


1338-1992 c 


18±4 


11+2 


0.29+0.05 


0.18+0.02 


(0.015,10%Z Q ) 


16.6 


12.2 


0.08 


0.06 



Table 4. Comparison between the Lya/Hell, Lyo/CIV, NV/Hell and NV /CIV measured (obs) and predicted (mod) values f or thr ee of 
the most extreme LAEs. Measurements have been taken from: "Humphrey 1 200 51 (Humphrey et al. 2006, in prep. ); ^Dey et al. ll997l: c De 
Brcuck et al. (in prep.). The best model characteristic parameters (U, Z) mo d f° r each object are shown in column 6. NV is observed to 
be much stronger than the models predict. 



once corrected for absorption), for which the two dimen- 
sio nal L ya spectrum appears very compact (De Breuck et 
al. 1 19991 ). or objects such that the aperture used to extract 
the spectra studied in this paper was optimized to isolate 
the nuclear emission. However, very high densities (>10 5 
cm -3 ) at large nuclear distances are implied by the enhanced 
Lya/Hell ratios measured in the extended gas of several ob- 
jects (see §4.1). 

An alternative possibility is that the gas that emits most 
of the line flux is distributed in very dense clumps, which 
can be at large di s tances from the nuclear region. As an ex- 
ample, iDev et all i|l997l ) find that 30% of the Lya flux in 
4C41.17 comes from extranuclear clumps in the scale of a 
few hundred pc. However, the implied photon ionizing lumi- 



nosities would be unrealistically high: Q ~10 



- 1 for 



r ~10 kpc, compared with ^severa l x 10 57 erg s _1 for z >2 
quasars (e.g. Heckman et al. Il991ah . 

We conclude that, in general, high densities do not pro- 
vide a natural explanation for the Lya excess in LAEs. 



5.3 Absorption and resonant scattering 

The simplest explanation for the larger Lya ratios in LAEs 
compared with objects with no Lya excess is that the line 
is less absorbed. This would also explain why Lya is also 
more luminous and has larger equivalent widths. It does 
not explain, however, why Lya/Hell is above the standard 
model predictions in some cases. 

A way to enhance Lya over other emission lines is by 
means of resonant scattering, with an absorbing medium 
distributed in very d ense and compact neutral dusty clumps 
(Hansen & Oh l2006l ). In such scenario, the Lya photons can 
be scattered by the gaseous surface of the neutral clumps 
without reaching the dust. All other emission lines and con- 
tinuum can penetrate the dusty cores, where they will suffer 
absorption. If the medium is quite clumpy, the Lya pho- 
tons will be able to escape unabsorbed. This can produce a 
noticeable increase of the line equivalent width and the ra- 
tio rel ative to other lines. IVillar-Martfn. Binette fc Fosburvl 
( 1996) proposed a similar explanation for the existence of 
pure Lya emitting regions in some r adio g alaxies. It has 
also been suggested by Hansen & Oh (2006!) to explain the 
intriguin g large EW Lya valu es measured in some Ly break 
galaxies. IVernet et all l|200ll) also found that in order to 
get an efficient grey scattering that explains the continuum 
polarization properties of HzRG, a high contrast clumpy 
medium is needed. 

In this scenario, the ambient medium around LAEs con- 
sists of dense neutral dusty clumps while objects with no 
Lya excess are surrounded by a more diffuse dusty medium 



where Lya is very efficiently quenched. However, as we ex- 
plain above (§5) LAEs are likely to be more strongly ab- 
sorbed than non-LAEs. 



5.4 Shocks 

So far we have investigated scenarios such that the same 
mechanism ionizes the gas in LAEs and non-LAEs, but the 
nebular properties are different. An alternative possibility 
is that the nebulae have similar properties, but the ionizing 
mechanism is different in both types of objects. 

Shocks could enhance the Lya emission via collisional 
excitation of the line, thanks to the heating effect and/or 
density enhancement. Two of the LAEs with the highest 
Lya excess (e.g. 1243+036 and 4C41.17) s how, in deed, clear 
sig ns of jet gas interactions (van Ojik et al. 1 19961 . Bicknell et 
al. 120001 ). LAEs in general tend to have small radio sizes and 
broad emission lines, also s uggestive o f jet-g as interactions 
(e.g. Humphrey et al. 120061 . Best et al. l200li ). 

The median values of the Lya luminosity for LAEs and 
non-LAEs arc 1.45 x 10 44 and 0.43 x 10 44 erg s" 1 respectively. 
Therefore, ignoring absorption, LAEs produce an excess of 
Lya luminosity given by L(Lya) e xc ~(1.45-0.43) x 10 44 erg 
s _1 =10 44 erg s _1 . There is, of course, a dispersion around 
L(Lya)exc, however, in order to simplify the argumentation 
that follows, we will consider the above value as the charac- 
teristic Lya excess for LAEs over typical radio galaxies with 
no Lya excess. 

If shocks are responsible for the Lya excess, it must be 
so that the other emission lines are not noticeably enhanced. 
Shock models with velocities in the range 200-1000 km s _1 , 
solar abundances and a gas precursor density of n =1 cm~ 3 
produce Lya/Hell in the range ~25-100 and Lya/CIV in the 
range ~10-100 for certain magnetic parameter values within 
the range B=10 _4 -10, which depend on the shock velocity 
(Mark Allen, private communication). For shock velocities 
larger than 400 km s _1 , Lya/H/3 is in the range 30-80. I.e., 
the shock cooling gas can be characterized by an emission 
line spectrum dominated by the Lya line. 

Shock models with Lya/Hell in the range 25-35 (as 
measured for the most extreme LAEs, i.e. those with the 
highest Lya excess) produce Lya/H/3 ~35. Therefore, the 
expected H/9 luminosity from the shocked gas L(H f3) shock is 
2.9xl0 42 erg s _1 . L(Hf3) s hock is related to the shock velocity 
in the clouds V w and the mass flow rate through the shock 
M by the following e quatio n (adapted from equation 4.4 of 
Dopita & Sutherland Il996h : 



L(HP) ah oek = 2.8xl0 3 



Vw 



M 



1 100 ferns" 1 J M Q yr 



— ergs 
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Figure 3. Density effects. Same diagnostic diagrams as in Fig.l with two model sequences for different U values. Solid line: U =0.05; 
dashed line: (7=0.7. The density n varies along each sequence between n=10 and 10 7 cm -3 or 10 6 cm -3 respectively. Density effects 
(coupled with a variation in U) can explain the location of all LAE in the diagnostic diagrams. However, unrealistically high densities 
(n ~10 5 cm" 3 ) are required to reproduce the line ratios of the highest Lya excess objects (see text). 



Assuming K„=500 km s -1 , this implies M <~10 4 M 
yr -1 . If the hot spots ad vance speed is in the range 
3xl0 3_4 kms _1 (e.g. Scheuer ll995h . the radio structures will 
need ~1.7xl0 6-7 yr to cross the ~50 kpc radius nebulae. 
The total amount of material consumed in the shock will be 
~1.7xl0 10-11 Mq, which is similar or larg er than the mass 
of the g iant nebulae (e.g. McCarthy 1 19931 . Villar-Martm et 
al. I2003T ) and is unrealistic. 

If we assume that L(H j3) s hock shows the same depen- 
dence on V w for values >500 km s _1 (although we do not 
know whether this is the case), extreme shock velocities as 
high as V w >5000 km s" 1 are required to obtain reasonable 
mass flow values <severalxl0 8 Mq. 

The continuum emitted by the shocked gas can also 
photoionize the pr ecurso r gas ahead of the s hock (e.g. Do- 
pita & Sutherland 1 19961 , Bicknell et al. l200fJ l. However, for 
photoionization of the precursor to enhance Lya efficiently 
over the other emission lines, very high densities and/or low 
metallicities would be required, as for standard photoion- 
ization (see above). As we discussed in §5.1 and §5.2, both 
scenarios have problems. 

Finally, we cannot explain either why Lya enhancement 
is not observed in the low z sample, where the extreme 
effects of sh ocks a re evident in ma ny ob jects (e.g. Villar- 
Martm et al. l2003l . Humphrey et al. 120061 '). 

We conclude that shocks do not provide a natural ex- 
planation for the Lya enhancement in the LAEs. 



5.5 Stellar photoionization 

An efficient way to enhance Lya over other emission lines is 
by means of stellar photoionization of the nebula, thanks to 
the soft ionizing continuum. 

Assuming that 67% of the ioniz ation processes end up 
in a Lya photon (Binette et al. 1 19931 ). L(Lya) exc ~10 44 erg 



s _1 corresponds to a total ionizing luminosity of Qlon ~10 55 
s" 1 absorbed by the nebula. 

A continuous burst of star formation with a star form- 
ing rate (SFR) of 200 M Q yr" 1 produces Q™° d =(0.9, 1.7, 
2.7) xlO 55 s _1 for ages of 1, 2 and 5 Myr respectively, con- 
sistent with Qion taking into account that a fraction of the 
photons can escape the nebula and Lya might have suf- 
fered some absorption Values up to a few thousands Mq 
yr -1 in HzRG have b een inferred from submm studies (e.g. 
Archibald et al. 120011 ). 

Let us compare the expected continuum flux at the red- 
shifted 1550 A with the observed values for two of the most 
extreme LAEs: 4C41.17 (z =3.79) and 1243+036 (z=3.57) 
for which / 15S o=l-3x 10" 18 and 3.4xl0" 19 erg s" 1 c m" 2 A 
s _1 respec tively in the observer's frame (Dey et al. Il997l ); 
Humphrey [20051 , Humphrey et al. 2006, in prep.). 

The expected monochromatic fluxes (nebular plus stel- 
lar, with the nebular component being <20% of the total 
continuum luminosity) of the bursts discussed above at the 
same wavelength are: (0.4, 0.8 and 1.6) xlO -18 erg s _1 cm -2 
A -1 for 1, 2 and 5 Myr respectively at z ~3.7, which are 
consistent with the measured values. 

For an instantaneous (rather than continuous) burst of 
star formation with the same IMF, a burst of age <2 Myr 
must have a mass 3.8xl0 8 Mq to reproduce Q abs ion (older 
ages require larger masses). The expected monochromatic 
flux at 1550 A (at z ~3.7) is /i 55 o=5.1xlO- 19 erg s" 1 cm" 2 
A _1 in reasonable good agreement with the continuum level 
of 1243+036. A mass ~10 9 Mq would be required to explain 
the continuum level of 4C41.17 (for this object, no contribu- 
tion from scattered continuum is expected, Dey et al. ll997l ). 

3 The escape fraction has very uncertain values. It can vary be- 
tween ~0.03 in some nearby starburst galaxies (Leitherer et al. 
1995) and > 0,5 in Lyman break galaxies (Steidel, Pettini &; Adel- 
bereer liuulT) . 
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It remains to be explained the trend for LAEs to show 
higher CIV/Hell compared with non-LAEs (median values 
1.7 and 1.0 respectively, see Table 3). This ratio is well above 
the standard power law model predictions for several objects 
(see Fig.l). 

A nebula photoionized by stars can be characterized 
by very large CIV/HellMOOO (Villar-Martfn et al. 2004). 
On the other hand, provided that the gas is highly ionized, 
CIV can have a non-negligible flux compared with Lya, with 
Lya/CIV>10 . Therefore, the UV emission line spectrum of 
such a nebula will be dominated by Lya and CIV, while 
the contribution to Hell (and NV) will be negligible. This 
is for instance the case of the Lynx arc, a gravitati onally 
lensed star forming galax y at z =3.36 (Fosbury et al. 120031 . 
Villar-Martin et al. l2004h . 

We show in Fig. 4 the UV line ratio predictions of mod- 
els which consider an increasing contribution of stellar pho- 
toionized gas relative to the AGN ionized gas (i.e. relative to 
the standard power law model predictions shown in Fig. 1). 
The uj parameter represents the ratio between the Lya flux 
emitted by the stellar ionized gas and that emitted by the 
AGN photoionized gas. Three values have been considered 
uj=Q (i.e. pure AGN photoionization, black dashed line), 0.7 
(red solid line) and 1.5 (green dotted line). For illustrative 
purposes, the stellar model predictions used here are those 
presented in Table 6 of Villar-Martm et al. 2004, correspond- 
ing to a single stellar age of 5 Myr and U —0.2 for the 
stellar photoionized gas. A Salpeter IMF in the mass range 
2-120 M© and an instantaneous star forming history were 
assumed. Lya/H/3 ~26 for this model. Fig. 4 shows that a 
varying contribution of stellar photoionized gas to the emis- 
sion line spectrum of LAEs is successful at reproducing the 
position of these objects in these three diagnostic diagrams. 

Stellar photoionization is, therefore, a promising expla- 
nation for the Lya excess in LAEs. A continuous star form- 
ing history with a star forming rate of ~200 M© yr _1 is 
very successful at explaining the excess of Lya luminosity 
and the observed level of continuum in at least the two most 
extreme LAEs. Since other objects show similar continuum 
levels and less extreme Lya excesses, it is likely that this 
assertion is valid for all LAEs. 

45% of LAEs require uj > 0: 1338-1941, 4C41.17, 
1243+036 and 0205+2242 (z >3) and 0303+3733 (z <3). 
The remaining LAEs are consistent with lu=0 within the er- 
rors (i.e. star formation is not required). However, since no 
information on the fraction of absorbed Lya flux is avail- 
able, except for 1755-6916 (which is not absorbed, Wilman 
et al. 2004), uj=0 is a lower limit. 



6 DISCUSSION 

There is additional evidence that supports stellar pho- 
toinization as the mechanism responsible for the Lya excess. 
Fig. 5 shows Lya/Hell (in log) vs. P(%), the polarization 
level of the optical continuum (UV rest frame) , for all radio 
galaxies at z >2 for w hich b oth have been measu red (Vernet 
et al. l200ll . Dey et al. Il997i . Jannuzi et al. Il995l . De Breuck 
et al. in prep.). There is a clear trend for objects with higher 
Lya/Hell values to show low polarization levels. In particu- 



lar, except for 1243+036, which breaks the trend, the LAEs 
(solid circles) show the lowest polarization levelQ 

The cor r elatio n in Fig. 5 is similar to that found by 
IVernet et al.l (|200ll ) on a sample of radio galaxies at 2< z <3 
between Lya/CIV and P(%) (all the objects in their sample 
are included in Fig. 5). Since the difference in polarization 
level between objects is most likely due to the diluting effect 
of young stars (Vernet et al. l200ll ). Fig. 5 strongly suggests 
that objects with higher Lya/Hell ratios contain a relatively 
more luminous stellar population. It is important to keep in 
mind that the range of Lya/Hell ratios must also be influ- 
enced by Lya absorption (it provides a natural explanation 
to why this ratio is well below standard PL model predic- 
tions in some objects). However, since some LAEs need an 
extra supply of soft ionizing photons to enhance the Lya 
emission, in spite of this uncertainty, Fig. 5 strongly sup- 
ports that stars provide such supply. 

Additional evidence for star formation has been found 
in several LAEs. This is the case of 4C 41.17 (Dey et al. Il997i . 
Dunlop et al Il994l . Chini & Kriigel Il994h and 1338-1942 
(Zirm et al. 120051 ). B3 0731+438, the only LAE in our study 
at 2 <3 for which optical polarization measurements exist, 
(Vernet et al. l200ll L shows null or very low polarization level 
(<2.4%), suggesting the presence of young stars. 

De Breuck et al. (in prep.) found a tentative trend (al- 
though more objects need to be studied to confirm this re- 
sult) for radio galaxies at z >3 to show lower optical contin- 
uum polarization level than radio galaxies in the 2< z <3 
range. This suggests that there is a relatively more lumi- 
nous young stellar population in the highest z radio galax- 
ies. Further evidence comes from the dramatic increase in 
the detection rate at submm wavelengths of z >2.5 radio 
galaxies (~75%) co mpared with z <2.5 radio galaxies (~15 
%) (Archibald et al. l200ll ). The authors also find that the 
average submm luminosity rises at a rate oc(l + z) 3 out to 
z ~4. Since most z >3 radio galaxies are LAEs, and pro- 
vided that the submm emi ssion is due to dust heate d by 
stars (e.g. Tadhunter et al. 120051 . Reuland et al. l2004h this 
adds further support to our interpretation that these contain 
more luminous young stellar populations. 

Enhanced Lya ratios, well above standard PL model 
predictions, have been measured at tens of kpc from the nu- 
cleus for the LAEs 4C41.17 and 1243+036 and some radio 
galaxies at 2< z <3 (§4.1), most of which are not LAEs. This 
suggests that star formation is spatially extended over large 
spatial scales. Spatially extended star formation in high z ra- 
dio galaxies has b een t entatively suggested by submm stud- 
ies (Stevens et al. 120031 '). For the low z objects, although the 
contribution of the stars to the Lya luminosity is not evi- 
dent in the spatially integrated emission line spectrum, this 
becomes apparent at tens of kpc from the nuclear region. 

Two more important issues need to be answered: Why 
is star formation more intense in LAEs? Was this class of 
objects much more common at higher (z Js3) redshifts and 
why? 



4 The polarization measurement for 1243+036 might not be reli- 
able, since it could be affected by strong sky residuals. The error 
bar represents a statistical estimate of the error based on noise 
properties and possible defects due to sky residuals are not in- 
cluded. 
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Figure 4. Effects on the UV line ratios due to the additional contribution of stellar photoionized gas. ui = ^y"* L represents the ratio 
between the Lya flux emitted by the stellar ionized gas and that emitted by AGN ionized gas. Same diagnostic diagrams as in Fig.l. The 
black dashed line represents the standard PL U sequence (see Fig.l ), i.e., the spectrum due to pure AGN photoionization (u>=0). The 
coloured lines represent models with an increasing contribution of emission lines from gas ionized by stars to the emission line spectrum, 
cj=0.7 and 1.5 for the red solid and green dotted lines respectively (see text for more details). U values for the PL sequence are shown 
for some models. The black solid arrows indicate the sense of increasing stellar contribution. A varying contribution of stellar ionized gas 
to the emission line spectrum of LAEs is very successful at reproducing the position of these objects in these three diagnostic diagrams. 
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Figure 5. log(Lya/HeII) vs. percentage of polarization of the op- 
tical continuum (UV rest frame) for all radio galaxies at z >2 for 
which both have been measured. Black solid symbols correspond 
to LAEs, while hollow symbols are objects with no Lya excess. 
Errorbars are shown when available. Notice the clear trend for ob- 
jects with the highest Lya/Hell ratios to show lower polarization, 
as expected if the Lya enhancement is due to stellar photoion- 
ization. The trend is broken by 1243+036, which shows too high 
P(%) for its large Lya/Hell ratio (but see text). 



The detailed spectroscopic study by iHumphrev et all 

( 2006) of a sample of 10 HzRG (z =2-3) show that small ra- 
dio sources exper ience stronger jet-gas interactions (see also 
Best et al. 1200 ll ) and contain more luminous young stel- 
lar populations. This could be the case of LAEs. In fact, 



this class of objects show in gen eral broad emiss i on line s 
and small radio sizes (se e Table 3; lDe Breuck et all (|2000bh . 
1)1501. De Breuck et al. 2000b) suggestive of jet-gas inter- 
actions. These are very clear in at least two o f the most 
extreme LAEs at z ^3: 4 C41.17 (Bicknell et al. |2000| ) and 
1243+036 (van Ojik et al. Il996l ) The star formation might 
have been induced by the interactions, although it is not 
clear that this process can form enough stars to dominate the 
rest-frame UV emission. A more attractive possibility is that 
the systems have recently undergone a merger event which 
has triggered both the star formation and the radio/ AGN 
activities. 

According to this study, we propose to use the Lya ra- 
tios (specially Lya/Hell) of high redshift radio galaxies and 
quasars to investigate the possible signature of star forma- 
tion. Although these ratios have been measured for many 
radio galaxies at different redshifts, this has been generally 
done for the spatially integrated spectra, where the effect is 
in general less clear. We propose that the high Lya ratios 
relative to other emission lines found in the extended emis- 
sion line regions of some h igh redshift radio-loud quasars 
(e.g. Heckman et al. Il991bl ) are a signature of star forma- 
tion, rather than low metallicities as often proposed. 

Star formation could also be responsible for the Lya 
emission detected in regions out of the reach of the quasar 
continuum and the radio structures i n som e high redshift 
radio galaxies (e.g. Villar-Martm et al. 2005, Reuland et al. 
l2003h . 

We argue next that the LAE fraction shows a genuine z 
evolution, not due to selection effects. Therefore, the radio 
galaxy phenomenon seems to have been more often associ- 
ated with a massive starburst at z >3 than at z <3. 
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6.1 Selection effects 

Our results suggest that LAEs were more abundant at higher 
redshifts. Whether this is a genuine z evolution depends on 
the real fraction of LAEs at different redshifts. 

An important problem we have when comparing the 
z <3 and z ^3 samples are the selection effects resulting 
from the search techniques used to find them and the tech- 
nical limitations of the spectrographs. 

The parent samples were selected based on radio prop- 
erties, meaning there is no a priori bias towards the opti- 
cal properties. However, higher power radio sources are ex - 
pected to emit stronger emission lines (Willott et al l 19991 ). 
so Malmquist bias on radio power could result on an optical 
bias towards strong line emitters. 

Is it possible that radio power is biasing our statistics? 
As De Breuck et al (2000) explain, the effects of Malmquist 
bias in their sample is alleviated for 2< z <4, with the addi- 
tion of several filtered surveys with lower flux density limits. 
We have compared the radio power at 1400 MHz (rest frame 
frequency, DBOOa) for the z ^3 and z <3 samples (Table 5). 
The median, average and standard deviation a are shown. 
According to the Kolmogorov-Smirnov test, the probability 
P (Table 5) that z >3 radio galaxies have different intrinisc 
radio power than the low z sample is very low (18%). So, 
we do not expect radio power selection criteria to introduce 
any bias. 

Could radio size be biasing our statistics? In samples to 
find HzRGs, the radio size has been used as an additional 
high redshift filter. As we discussed above, since small radio 
sources tend to have more intense star formation, this selec- 
tion criteria could introduce some additional bias. However, 
the 2< z <3 and z ^3 samples are not different in terms of 
their radio sizes (see §4, Table 2) . On one hand this selection 
criteria was applied for only ~10% of the sources in Table 1. 
More importantly, 0140+3253 is the only z ^3 radio galaxy 
from such samples. The other 12 z Js3 radio galaxies have 
been found without radio size constraints. This selection can 
therefore not explain the different LAE fractions at z <3 and 
z^3. 

We conclude that no bias is introduced due to radio 
selection criteria. 

Our work is affected by the need to detect strong emis- 
sion lines since a) this was required in the original sample 
to measure the redshift b) we need at least Lya and Hell 
to be observed to determine whether an object shows Lya 
excess. What is clear is that, due to their Lya excess com- 
pared with non-LAEs, we are not likely to have missed low z 
LAEs, taking also into account that at low z, CIV and Hell 
are in a sensitive part of the spectrographs. There are 26 
objects in De Breuck et al (DBOOa DB01) at 2< z <3 which 
could not be classified as LAEs or non LAEs (i.e., Lya/Hell 
not available). An extreme, unlikely case of all these objects 
being LAEs would be required for the fractions of LAEs in 
the two z samples to become similar (46% vs. 54%) On the 
contrary, we would rather expect to have missed non-LAEs 
in the low z sample with weak emission lines. This would 
make the fraction of LAEs in the low z sample even lower. 

Is it possible that the sources with no CIV and Hell 
are systematically biasing our statistics in the z ^3 sample? 
In the De Breuck et al. (2000) sample, there are 22 objects 
at z ^ 3; 11 of these could not be classified as LAE's or 
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Table 5. Comparison between the rest frame radio power at 1400 
MHz (in units of 10 35 erg s" 1 Hz" 1 ) for the 2< z <3 and z ^3 
samples. Radio power value taken from DBOOa. The median, aver- 
age and standard deviation values are shown. P is the probability 
that the two samples are drawn from different parent populations 
according to the Kolmogorov-Smirnov test. This tests implies that 
there is no significant difference in intrinsic radio power between 
both samples. 

non-LAE's. Only 6 of these had Hell outside the observed 
wavelength range. The remaining 5 sources had too shallow 
limits on the Lya/CIV or Lya/Hell ratio to classify them. 
Even in the extreme case that these 11 objects were non- 
LAEs, the fraction of LAEs at z ^3 would be 7 out of 22, 
or 32%, still significantly larger than in the low z sample. 

Is it possible that non-LAEs are systematically missing 
from samples of HzRGs? In the sample of ultra steep spec- 
trum selected HzRG candidates of De Breuck et al (2001), 
13 (35%) sources did not yield a redshift. From these, six did 
not show any emission, while the other seven had only con- 
tinuum emission down to the blue end of the spectrograph, 
but no emission lines. The latter sources have to be at z < 2 
(and therefore of no interest in the present work) based on 
the absence of a continuum discontinuity across Lya. 

The former sources could be at z ^2. The emission lines 
may be in the observed wavelength range, but too faint to 
be detected. If they were objects with Lya in the observed 
wavelength range, it would have been more likely that they 
would be non-LAEs, due to their intrinsically fainter Lya 
emission. In any case, these sources would add at most six 
objects to our exisiting HzRG sample. In the extreme, un- 
likely case that these were z ^3 non-LAEs, the fraction of 
high z LAEs would change from 54% to 37%, still signif- 
icantly larger than in the low z sample. If these, together 
with the 11 2 ^3 non-classified objects were z ^3 non-LAEs, 
the fraction would be 25% vs. less than 10% expected in the 
low z sample. 

We conclude that although the fraction of LAEs may 
be incompletely determined, both at 2< z < 3 and at z ^ 3, 
the large difference between both fractions is likely to be 
real and not a consequence of selection effects. An extreme, 
unlikely situation would be required to change significantly 
the large difference between the LAE fractions: most non- 
classified objects at 2< z < 3 should be LAEs and most 
non-classified objects at z ^3 should be non-LAEs. 

In addition, we note that all z > 3 LAE's are at z > 3.5. 
If we compare the z < 3.5 and z > 3.5 sources, the fraction of 
LAE's are 4 out of 52 (8%) and 7 out of 9 (78%), respectively. 
It thus seems impossible to explain such large differences as 
due to selection effects. 



7 SUMMARY AND CONCLUSIONS 

The behaviour of the Lya/Hell, Lya/CIV ratios and the 
Lya line luminosity has been investigated in a sample of 48 
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radio galaxies at 2< z <3 and 13 radio galaxies at z ^3 for 
which Lya/Hell and/or Lya/CIV measurements are avail- 
able. 

As found bv lDB OOa. z Js3 radio galaxies tend to show 
larger Lya ratios (2.2 and 1.5 times higher for Lya/Hell and 
Lya/CIV respectively) and higher (3 times) Lya luminosi- 
ties than the low z sample. The Lya ratios are consistent 
with or above case B predictions of standard power law pho- 
toionization models in at least 54% (7 out of 13) of the z ^3 
objects. This is the case for only 8% (4 out 48) of radio galax- 
ies at 2^ z <3. We refer to all these objects with unusually 
enhanced Lya emission as Lya excess objects or LAEs. 

LAEs show Lya/Hell, Lya/CIV and Lya luminosity 
values which are ~2-3.5 times higher than non-LAEs. They 
have radio sizes ~4.6 times smaller. 36% of LAEs show 
Lya/Hell above the standard model predictions. Since in 
general no absorption has been taken into account, this frac- 
tion could be larger. 

Several possibilities have been investigated to explain 
the Lya excess in LAEs: stellar photoionization, shocks, low 
metallicities, high densities and absorption/resonant scat- 
tering effects. All scenarios but the first can be rejected 
with confidence: we propose that a population of young stars 
generate the extra supply of soft ionizing photons necessary 
to produce the observed Lya excess in LAEs. Star form- 
ing rates ~200 Mq yr _1 are required. The enhanced and 
intense star formation activity in LAEs could be a conse- 
quence of a recent merger event which has triggered both 
the star formation and the radio/ AGN activities. This in- 
terpretation is strongly supported by the clear trend for ob- 
jects with lower optical continuum polarization level to show 
larger Lya/Hell ratios. It is further supported by the ten- 
tative trend found by other authors for z >3 radio galaxies 
to show lower UV-rest frame polarization levels, or the dra- 
matic increase on the detection rate at submm wavelengths 
of z >2.5 radio galaxies. 

The models imply that star formation is occuring in 
45% (5 out of 11) of LAEs. Since absorption has not been 
taken into account, this is a lower limit. 

The finding of enhanced Lya emission at tens of kpc 
from the nuclear region in some objects suggests that star 
formation can be extended over spatial scales of tens of kpc. 

We argue that, although the fraction of LAEs may be 
incompletely determined, both at 2< z < 3 and at z ^ 3, 
the much larger fraction of LAEs found at z ^3 is a genuine 
redshift evolution and not due to selection effects. There- 
fore, our study suggests that the radio galaxy phenomenon 
was more often associated with a massive starburst at z >3 
than at z <3. In other words, powerful radio galaxies (and, 
according to the unification model, powerful radio quasars), 
appear in more actively star forming galaxies at z >3 than 
at z <3. 

Enhanced Lya ratios in the extended gas of some z <3 
radio galaxies suggest that star formation is also present in 
these objects, although at lower levels. Star formation could 
also be responsible for the Lya emission from regions out of 
the reach of the quasar continuum and the radio structures 
detected in some high redshift radio galaxies. 

We propose to use the Lya ratios (specially Lya/Hell) 
of high redshift (z >2) radio galaxies and the extended gas in 
radio loud and radio quiet quasars to investigate the possible 
signature of star formation. 
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